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ABSTRACT 

We investigate the uncertainty in surface abundances and yields of asymptotic giant 
branch (AGB) stars. We apply three different mass loss laws to a 1.5 star of 
metallicity Z = 0.008 at the beginning of the thermally pulsing asymptotic giant 
branch (TP-AGB) phase. Efficient third dredge-up is found even at very low envelope 
mass, contrary to previous simulations with other evolution codes. We find that the 
yield of carbon is uncertain by about 15% and for most other light elements the 
yield is uncertain at the level of 20-80%. For iron group elements the uncertainty 
varies from around 30% for the more abundant species to over a factor of two for the 
less abundant radioactive species, like ^"Fe. The post-AGB surface abundances for 
this mass and metallicity are much more uncertain due to the dilution of dredged-up 
material in differing envelope masses in the later stages of the models. Our results are 
compared to known planetary nebula (FN) and post-AGB abundances. We find that 
the models are mostly consistent with observations but we are unable to reproduce 
observations of some of the isotopes. 
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1 INTRODUCTION 

Stars of between about 1 and 8 M© experience a phase of un- 
stable double shell burning known as thermal pulses whilst 
on the asymptotic giant branch (AGB). During this phase, 
the star experiences significant mass loss, which strips away 
the envelope to leave behind a degenerate core of carbon 
and oxygen. The mechanism for this mass loss is unknown 
and this lack of knowledge is a major source of uncertainty 
in calculations of thermally pulsing asymptotic giant branch 
(TP-AGB) evolution. Mass loss allows the nucleosynthetic 
products of AGB evolution to be returned to the inter- 
stellar medium and, because low-mass stars are formed far 
more abundantly than those of higher mass (|Kroupa et al.l 
1 19931 ). allows AGB stars to contribute significantly to galac- 
tic chemical evolution. Differences in mass-loss history may 
affect the evolution and nucleosynthesis of a star and hence 
lead to different yields. 

Mass loss also has important effects on the evolution of 
AGB stars. In the lowest-mass AGB stars it is possible for 
mass loss to strip the star of its envel ope before third dred ge- 
up (TDUP) has a chance to occur (|Karakas et al. |[2003). If 
TDUP does occur its efficiency can still be affected by mass 
loss. In intermediate-mass AGB stars the occurrence of mass 
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loss can inhibit the occurrence of hot-bottom burning in the 
later stages of evolution, while still allowing TD UP to occur. 
This allows the formation of bright carbon stars (|Frost et al.l 
Il998l ). 

Of the early calculations of A GB evolution, few used 
any mass loss. Those that di d (e.g. Schonberneij[l979h em- 
ployed the Reimers' formula (|Reimer sl l 19751 ) 



r p 

M = -4 X 10~"»7— Moyr" 



(1) 



where L, R and M are the star's luminosity, radius and 
mass in solar units. The value of the free parameter 77 is 
usually t aken as .4 for red giants. However, it was recog- 
nized by iRenzinil (|l98lh that the mean mass-loss rate re- 
quired to produce a typical planetary nebula was around 
3 X 10~^ MQyr~^ . He coined the term superwind as this value 
was significantly greater than those given by the Reimers 
formula. Many calculations still use the Reimers formula, 
generally with 0. 4 < 7j < 3, although t? = 10 has been used 
in extreme cases (jStraniero et al ] |l997l . l2003l ). 

In an attempt to produce a more physically consistent 
picture. I Vassiliadis fc Wood (1993) produced a mass- loss re- 
lation based on observations of OH/IR stars and other pul- 
sating, dust-enshrouded AGB stars. They linked the mass- 
loss rate, M in M0yr~^, of an object to its Mira pulsational 
period, P (in days), via the relations: 
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log M ■ 
and 



-11.4 + 0.0123P 



M 



(2) 



(3) 



where «exp ~ —13.5 + 0.056P is the expansion velocity of 
the wind far from the star in kms~^ and c is the speed of 
light. The mass loss is taken to be the smaller of the two 
values calculated with the above expressions. This law has 
the advantage of providing for a superwind phase which can 
remove a star's envelope over the course of a few pulses. 
This mass loss law h as been used in recent computations by 
iKarakas etaP (|2002l l. 

One other mass-loss relation has al so been used in ex- 
tant calculations, that of iBlockerl l|l995h . It is grounded in 
theory, rather than empirically determined, based as it is 
on simulations of shock-driv en winds in the atmospheres of 
Mira-like stars (|Bowenlll988l ). It was developed to reflect the 
strong increase of mass loss during the A GB and to be ap- 
plicable to stellar evolution calculations l|Blockeilll995l ). It 
is similar to the Reimers' prescription but with a steeper 
dependence on luminosity and a shallower dependence on 
mass. The formula is expressed as: 



M = 4.83 X 10" 



M 



Mr 



(4) 



Because of its use of the Reimers' formula, this law also 
suffers from the problem of having a free parameter that 
must be determined. Recent calculations using Blo cker's re- 
lation include those of iHerwig fc AustinI l|2004h and lVentural 
ll2004h . 

In this work, we address the uncertainty in the evolution 
and nucleosynthesis of a low-mass TP- AGB star due to our 
lack of understanding of the mass loss mechanism. 



2 THE STELLAR EVOLUTION CODE stars 

We use the stellar evolution code STARS which is a vari- 
ant of the evolution code originally devel oped by [ Eggletonl 
and updated by many authors (e.g. lPols et al.lll995l ). 
The ve rsion used here includes the A GB-specific modifica- 
tions of lStancliffe. Tout fc Po5 ll2004h, together with the up- 
dated opacity tables of'Eld ridge fc Tou^ (|2004l ') . In order to 
study light element nucleosynt hesis in detail we emp loy the 
nucleosynthesis subroutines of lStancliffe et al] (|2005l ) which 
cover isotopes from deuterium to sulphur and important iron 
group elements - a total of 40 isotopes. 

A 1.5 M0 model of metallicity Z = 0.008 was evolved 
from the pre-main sequence using 199 mesh points. Initial 
abundances are solar-scal ed and taken from the work of 
lAnders fc Grevessel l| 19891 ). On the first giant br anch and 
on th e early asymptotic giant branch (E-AGB) iReimerj 
1 19751 ) mass loss was used with 77 = 0.4. At the end of 
core helium burning, the model was remeshed to 499 mesh 
points to facilitate the transition to the high resolution 
(999 mesh points), A GB-specific mesh spacing function of 
IStancliffe et al. I l|2004h . Convective overshooting was not em- 
ployed during any stage in the evolution. 

On the TP- AGB, we use 77 = 1 when applying Reimers' 
mass los s. When applying Blocker's law, we us e rj — 0.02, 
following I Ventura. D'Antona fc MazzitellH (|2000l ) who chose 
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Figure 1. Plot of the total stellar mass as a function of time from 
the beginning of the evolution run. The model using Reimers' 
mass loss is the solid line. The dashed line is the model with 
Blocker's mass loss and the dotted line is the model with Vassil- 
iadis and Wood (VW) mass loss. 



this value based on calibration to the lithium luminosity 
function in the Large Magellanic Cloud. 



3 RESULTS 

The three mass-loss laws lead to different mass loss histories 
for each of the models. The evolution of stellar mass with 
time is shown in Figure [1] The model with Reimers' mass 
loss has appreciable mass loss from the beginning of the run. 
The rate of mass loss slowly increases with time before the 
envelope is completely removed. The model with Blocker's 
mass loss only begins to show appreciable mass loss after 
about 5 X 10^ yr. The mass loss then increases more rapidly 
as the model evolves. The model with Vassiliadis and Wood 
(hereafter VW) mass loss only shows appreciable mass loss 
once the superwind phase begins at around lO'' yr. 

When the model using VW mass loss reaches an enve- 
lope mass of around 3 x 10"'^ M0, the star experiences its 
last thermal pulse. Because the envelope mass is very low, 
the pulse occurs when the star is slightly further to the blue 
in the HR diagram (see Figure [^. Th is is an AGB final 
thermal pulse (AFTP, see lBlockeij|200ll ). 

Details of the evolution of the three models are pre- 
sented in Tables [1] [2] and |3] The Reimers' model undergoes 
the most thermal pulses but all models have just six ther- 
mal pulses followed by episodes of third dredge-up. Dredge- 
up begins slightly later in the the Reimers model, when the 
core mass has reached 0.602 Mq compared with 0.595 Mq 
and 0.594 Mq in the Blocker and VW models respectively. In 
terms of the total amount of intershell material dredged-up, 
all the models are comparable. The Reimers' model dredges 
up 0.0179 Mq of material while for the Blocker and VW 
models, 0.0171 Mq and 0.0176 Mq of intershell material are 
dredged up respectively. Because the Reimers model typi- 
cally has a lower envelope mass for a given thermal pulse, 
when dredge-up occurs the material that is carried to the 
envelope does not get as diluted as it does in the other two 
models so a higher C/O ratio tends to be obtained (exclud- 
ing the AFTP as this is a slightly different phenomenon). 
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TP M, Afenv Mh Tip log(L™/L0) AA/h AAfnup A C/O 

(Mo) (Mq) (Mq) (104 yr) (M0) (Mq) 



1 


1 


36357 


0.79285 


0.57072 




4.63556 





01112 









324 


2 


1 


30354 


0.71765 


0.58183 


7.94 


6.42507 





00405 









324 


3 


1 


27895 


0.68858 


0.58587 


9.57 


6.48659 





00449 









324 


4 


1 


24563 


0.64962 


0.59035 


9.67 


6.89601 





00566 









324 


5 


1 


20506 


0.60288 


0.59599 


9.59 


7.04725 





00619 









324 


6 


1 


15642 


0.54764 


0.60203 


9.30 


7.20978 





00675 


0.00134 


0.199 





553 


7 


1 


09658 


0.48161 


0.60744 


9.02 


7.40975 





00753 


0.00252 


0.335 


1 


086 


8 


1 


02219 


0.40142 


0.61245 


8.80 


7.60488 





00832 


0.00373 


0.448 


1 


993 


9 





92994 


0.30389 


0.61704 


8.52 


7.79146 





00902 


0.00303 


0.336 


2 


860 


10 





82767 


0.19597 


0.62303 


7.84 


7.81919 





00867 


0.00308 


0.355 


4 


148 


11 





71431 


0.07700 


0.62862 


7.21 


7.88218 





00872 


0.00422 


0.484 


8 


072 



Table 1. Details of the model computed using Reimers' mass loss with r] = 1. The data are TP — the thermal pulse number, Af, - the 
total stellar mass, Afenv - the envelope mass, A^h — the hydrogen free core mass. Tip - the interpulse period, LjJ^'' - the peak luminosity 
from helium burning, AA/jj - the hydrogen free core mass growth during the interpulse, AAfouP ^ th^ mass of material dredged up, 
A = AAfDup/AAfjj - the dredge- up efficiency and C/O - the surface carbon-to-oxygen ratio by number. 



TP 


M, 


Afcnv 


Mh 


Tip 


log (Lg--/Lo) 


AA/h 


AAfnup 


A 


C/O 




(Mo) 


(Mq) 


(Mo) 


(104 yr) 




(Mo) 


(Mo) 






1 


1.36357 


0.79285 


0.57072 




6.48915 


0.01468 






0.324 


2 


1.34865 


0.75917 


0.58538 


9.68 


6.31532 


0.00410 






0.324 


3 


1.33271 


0.73750 


0.58947 


9.82 


6.91354 


0.00575 






0.324 


4 


1.30801 


0.70657 


0.59518 


9.88 


7.05616 


0.00626 


0.00054 


0.086 


0.362 


5 


1.27055 


0.66274 


0.60090 


9.48 


7.23129 


0.00691 


0.00162 


0.234 


0.589 


6 


1.21258 


0.59870 


0.60619 


9.21 


7.43439 


0.00769 


0.00274 


0.356 


1.048 


7 


1.12323 


0.50370 


0.61114 


8.93 


7.61591 


0.00839 


0.00367 


0.437 


1.754 


8 


0.98724 


0.36246 


0.61586 


8.59 


7.75249 


0.00892 


0.00451 


0.506 


2.873 


9 


0.77352 


0.14381 


0.62027 


8.26 


7.91912 


0.00945 


0.00402 


0.425 


5.104 



Table 2. Details of the model computed using Blocker's mass loss with r) = 0.02. The columns are the same as in Table [T] 

Recent work bv lStancliffel | 2006t ) has shown that this is not 
due to the use of a simultaneous method of solution of the 
equations of stellar structure and evolution. 

The physical cause for the occurrence of TDUP has re- 
ceived li ttle a ttention in recent years. It is briefly touched 
upon bv llbeni (jlOTCj) who states: 

The entropy content of the convective envelope above the 
hydrogen-helium interface is immense relative to the entropy con- 
tent of the helium-rich region below the interface. Hence, the ab- 
sorption of comparable amounts of energy by each region corre- 
sponds to a much, much lower increase in entropy in the convec- 
tive region above the interface than in the radiative region below 
the interface. The only way in which the continuity of the en- 
tropy parameter S can be maintained across the interface is for 
the interface to move inward in mass. 

This description is rather brief. 

We offer the following description of why TDUP occurs. 
Energy from helium burning causes expansion throughout 
the star. This expands the hydrogen burning shell and ex- 
tinguishes it. The outermost envelope initially contracts be- 
cause it is no longer being supported by H-burning - this 
explains the immediate dip in both surface luminosity and 
radius seen after a thermal pulse. Eventually, energy from 
He-burning filters through to the these regions and expan- 
sion resumes. 

This expansion has the following effects. The expan- 
sion causes a drop in temperature and hence the opacity 



4.5 




5 4.5 4 3.5 

loglO Teff/K 

Figure 2. Ifertzsprung-Russell diagram for each of the three 
models. The model with VW mass loss experiences an AGB final 
thermal pulse. 

3.1 Dredge- up at low envelope mass 

One important feature of these models is the occurrence 
of third dredge-up even at very low envelope mass. Pre- 
vious studies have suggested that TDUP s hould cease be- 
low envelope masses of around 0.5 M© (IStraniero et al.l 
I2OO3I ). Calculations with the STARS code have given deeper 
dredg e-up in low-mass star s compared with other simula- 
tions (IStancliffe et al.ll2005l ). The reason for this is unclear. 
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TP Afenv Mh Tip log(L™/L0) AMh AA/dup A C/O 

(Mo) (Mq) (Mq) (104 yr) (M0) (Mq) 



1 


1.36357 


0.79285 


0.57072 




6.45469 


0.01467 






0.324 


2 


1.36337 


0.77396 


0.58537 


9.45 


6.30073 


0.00403 






0.324 


3 


1.36293 


0.76796 


0.58939 


9.59 


6.86964 


0.00557 






0.324 


4 


1.36047 


0.75949 


0.59494 


9.64 


7.00581 


0.00603 


0.00055 


0.091 


0.359 


5 


1.35569 


0.74849 


0.60042 


9.33 


7.18123 


0.00678 


0.00175 


0.258 


0.587 


6 


1.34349 


0.73042 


0.60545 


9.10 


7.39886 


0.00762 


0.00314 


0.412 


1.047 


7 


1.28632 


0.66795 


0.60993 


8.91 


7.60311 


0.00844 


0.00393 


0.466 


1.652 


8 


0.98573 


0.36262 


0.61444 


8.43 


7.68659 


0.00867 


0.00432 


0.498 


2.798 


9 


0.63154 


0.00362 


0.61879 


8.22 


7.83213 


0.00914 


0.00403 


0.441 


28.01 



Table 3. Details of the model computed using Vassiliadis and Wood mass loss. The columns are the same as in Table [T] Note that the 
last pulse is an AGB final thermal pulse. 



Grad J 
Graded 









0.62786 0.62787 0.62788 0.62789 0.6279 



Figure 3. The adiabatic and radiative temperature gradients 
as a function of mass at the base of the convective envelope. 
The adiabatic gradient, Vadi is the solid line. The dashed line 
represents the radiative gradient, Vrad- This profile is taken from 
close to the beginning of TDUP in the last thermal pulse in the 
VW model. The plot is not smooth due to the number of decimal 
places that the output file uses to print the mass co-ordinate. 



in the envelope increases. This causes a sharp rise in the 
radiative temperature gradient, Vrad, pushing it above the 
relatively constant adiabatic gradient, Vad- This makes the 
region convective and leads to the occurrence of TDUP. The 
increase in entropy described by Iben is a consequence of the 
expansion, not a cause of third dredge-up. 

Should we expect TDUP to continue even at low enve- 
lope masses? In the picture outlined above, the answer would 
seem to be 'yes'. If sufficient energy from helium burning is 
dumped into the envelope, then TDUP should occur. 

Why do the models presented here give such efficient 
dredge-up at low envelope mass? Most importantly, the ra- 
diative temperature gradient does not show a discontinuity 
as it approaches the adiabatic gradient from the convective 
side of the border, as can be seen in Figure |3l This is an im- 
portant thing to achieve in calculations of third dredge-up. 
A discontinuity in the radiative gradient would correspond 
to an unstable Schwarzchild boundary: mixing of material 
across this boundary would result in the moveme nt of the 
locat ion of this boundary (see the discussion in iMowlavil 
I1999I . for example). It is thus necessary to ensure that a sta- 
ble boundary is obtained, as is done here, to ensure that 
dredge-up is accurately calculated. 



It may be suggested that the depth of TDUP depends 
on the choice of the mixing coefficient at the edge of a con- 
vective zone. The STARS code employs an arithmetic mean to 
compute t he mixing between tw o mesh points (see the dis- 
cussion m IStancliffe et"al]l2004 ). It can be argued that this 
choice amounts to a form of convective overshooting as we 
are mixing material from beyond the formal (Schwarzchild) 
boundary for convective stability. However, we are mixing 
only one zone beyond the Schwarzchild boundary and this 
is much smaller that the scale of convective overshoo ting 
employed in a code like that of IStraniero etHI (j2003l ). as 
the mass resolution around the base of the convective is of 
the or der of 10~® Mq. We also point out that lCristallo et al.l 
(|2004 ) did not find efficient dredge-up at low envelope mass 
despite the inclusion of convective overshooting. The reason 
why we find deep dredge-up at low envelope mass without 
the use of convective overshooting while others do not (even 
though they include overshooting) is not understood. 

3.2 Nucleosynthesis 

Before discussing the nucleosynthesis of these models it is 
necessary to point out a deficiency. Because there is no con- 
vective overshooting (or other extra mixing process such as 
rotationally-driven mixing or internal gravity waves) used 
in the models, no carbon-13 pocket is formed. It is the re- 
action ^^C(Q:,n)^^0 that is thought to be responsible for 
the s-process in low-mass AGB stars. This should be borne 
in mind when reactions involving neutron captures are dis- 
cussed below. 

Light element yields and final surface abundances for 
the three models are presented in Tables 111 and [^. For ease 
of comparison, they are also displayed in Figures |4] and [5] 
The signatures of nucleosynthesis along the whole TP- AGB 
are displayed most clearly in the model with Reimers' mass 
loss. This is because it typically has a lower envelope mass 
for each thermal pulse and so the dredged-up material suffers 
less dilution in the envelope. This explains why the surface 
abundances of some isotopes are significantly higher in this 
model compared with the model with Blocker mass loss. 
Because there is almost no envelope remaining when the VW 
model experiences its AFTP, its final envelope composition 

^ Yields and abundances for all the isotopes in the nucleosynthe- 
sis network can be found in the Appendix. 
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Figure 4. Net yield in solar masses of selected isotopes from our models. Squares are yields from the model with Reimers' mass loss, 
circles are from the m odel with Block er's mass loss and triangles are from the model with VW mass loss. For comparison, we have also 
included the yields of iKarakasI | |2003|) . which are displayed with crosses. Only those isotopes with a positive net yield are displayed - 
hence the absence of yields for the isotopes ^^C, ^^O and ^^Mg from the Karakas model: her net yields for these isotopes are all negative. 



is dominated by the nucleosynthesis signatures of this last 
pulse cycle. 

Figure |4] also plots the yields of a 1.5 star o f metal- 
licity Z = 0.008, as calculated by iKarakaa l|2003h . along 
side our yields. The Karakas model was calculat ed using 
the M ount Stromlo Stellar Structure Program fsee iKarakasI 
|2003| . and references therein). It should be noted that the 
K arakas model uses the init ial ^^C, ^''N and abundances 
of I Russell fc Dopital l|l992l ). rather than solar-scaled abun- 
dances as we have used. For isotopes like ^'*N, whose abun- 
dance is mostly determined by evolution prior to the TP- 
AGB, we obtain very similar yields to those of the Karakas 
model. However, there are significant differences (of sev- 
eral orders of magnitude) between the yields of the Karakas 
model and the yields of our three models for those isotopes 
that are produced in the intershell of an AGB star. Most 
notable are the yields for ^^C, ^"^O and ^'^Mg. For these iso- 
topes we obtain a positive net yield (see equation[5]below for 
a definition), whereas her model has a negative net yield. We 
attribute these differences to the fact that all of our models 
show quite efficient third dredge-up whereas hers does not. 

When the models first reach logjp Toh = 5.0, the sur- 
face carbon abundance (by mass) of the model with Reimers' 
mass loss is 2.636 x 10~^ compared to 1.602 x 10~^ in the 
model with Blocker's mass loss and 1.579 x 10~^ in the model 
with VW mass loss. However, because all three models un- 
dergo comparable amounts of dredge-up and lose roughly 



the same total mass, their carbon yields (i.e. material re- 
turned to the interstellar medium) are comparable, with a 
spread of just 6 x 10~^ M©. This amounts to an uncertainty 
of about 15% associated with the mass-loss law used. All 
yields described herein are net yields, i.e. 

p = [ M(t)[X4t)-Xs{0)]dt, (5) 
Jo 

where M{t) is the mass loss rate at time t, Xa{t) is the 
mass fraction of species s at time t and Xa{0) is the initial 
mass fraction of species s in the material from which the star 
formed. 

The post- AGB abundance shows an uncertainty 
of around 67%. Oxygen is produced in the intershell dur- 
ing helium burning and brought to the surface during third 
dredge-up. It is only a minor component of the intershell 
(at around 2% by mass) and is little affected by reactions 
in the H-burning shell. Hence there is little variation in the 
surface abundance (see Table except in the VW model 
when the last pulse makes a significant contribution to the 
oxygen content of the envelope. The yield is well de- 
termined, with an uncertainty of 6.6%. This is because the 
models all dredge-up a similar quantity of He-burning ashes. 

The pathway to ^^F production in AGB s tars is a com- 
plicated one. It is produced from ^''N via l|Gorielv et al.l 
Il989ll 
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Figure 5. Plot of the logarithm of the ratio of the final surface abundance to the initial surface abundance for selected isotopes in the 
nucleosynthesis network. Squares are from the model with Reimers' mass loss, circles are from the model with Blocker's mass loss and 
triangles are from the model with VW mass loss. 



Mass loss 




16o 


19 p 


21 Ne 


22Ne 


23Na 


2BMg 


26Mg 


26A19 


Reimers 
Blocker 
VW 


4.525(-3) 
4.449(-3) 
3.884(-3) 


1.197(-4) 
1.155(-4) 
1.234(-4) 


3.730(-7) 
3.428(-7) 
2.367(-7) 


3.571(-8) 
3.091(-8) 
1.522(-8) 


2.603(-4) 
2.266(-4) 
1.646(-4) 


6.301(-6) 
5.419(-6) 
4.016(-6) 


2.1.38(-6) 
1.506(-6) 
9.254(-7) 


1.255 (-6) 
9.753(-7) 
6.970(-7) 


6.104(-8) 
6.081(-8) 
5.498(-8) 


Difference 
(%) 


6.41(-4) 
15.0 


7.90{-6) 
6.6 


1.363 (-7) 
42.9 


2.05(-8) 
75.1 


9.57(-5) 
44.1 


2.29{-6) 
43.5 


1.22{-6) 
79.6 


5.58(-7) 
57.2 


6.06{-9) 
10.3 



Table 4. Net yields (in solar masses) of selected light isotopes. All yields are expressed in the form n(m) = n X 10"*. The percentage 
difference is calculated using the mean of the three yields as a reference point. 



Mass loss 


12c 


16o 


19p 


2lNe 


22Ne 


23Na 


2BMg 


26Mg 


26A1S 


Reimers 


2.636(-2) 


4.393(-3) 


2.859(-6) 


1.947(-6) 


1.884(-3) 


5.112(-5) 


4.514(-5) 


4.017(-5) 


2.031(-7) 


Blocker 


1.602(-2) 


4.179(-3) 


1.418(-6) 


1.777(-6) 


9.074(-4) 


2.955(-5) 


3.384(-5) 


3.471(-5) 


1.314(-7) 


VW" 


l..->7<)(-l) 


7.,S()7(-:-!) 


l.i3()(-.-)) 


3.12S(-(i) 


l).')8.^)(-:-!) 


1.772(-1) 


1.20()(-i) 


7.70:-!(-.-)) 


i).(i71(-7) 


Difference 


1.42(-1) 


3.69(-3) 


1.29(-5) 


1.35(-6) 


9.08(-3) 


1.47(-4) 


8.61(-5) 


4.23(-5) 


8.36(-7) 


(%) 


212 


67.3 


208 


59.1 


719 


171 


130 


83.6 


193 



Table 5. Surface abundances by mass fraction for selected light isotopes when the models reach logj^Q T(,ff = 5.0 for the first time. All 
abundances are expressed in the form n{m) = n x 10'". The percentage difference is calculated using the mean of the three abundances 
as a reference point. 
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"N(a, 7)^«F(/3+)^«0(p, a)^^N(a, -y^^F. (6) 

In order to get p-capture on to in the H-depleted inter- 
shell, the ^^C(a,n)^®0 reaction must be active in order to 
produce neutrons so that the reaction ^'*N(n, p)^''C can oc- 
cur. Uncertainties in fluorine production in AGB stars due 
to uncertainties i n the reaction rates involved were studied 
bv lLugaro et al. They found an uncertainty of about 

50% in the yields from their low mass models. This is compa- 
rable to the uncertainty found here (43%) due to the uncer- 
tainty in mass-loss rate. The post- AGB surface abundance 
of ^^F is found to be highly uncertain (by a factor of 3), with 
the model with VW mass loss having a much higher post- 
AGB surface abundance. The model with Reimers' mass loss 
also has a higher surface abundance than the Blocker's mass 
loss model. This is because in the VW model sequence, the 
dredged-up material suffers less dilution in the low mass of 
envelope remaining at the time of the last thermal pulse. 
The sequence using Blocker's mass loss has a more massive 
envelope at its last thermal pulse, so the dredged-up mate- 
rial suffers more dilution. 

Of the Ne-Na cycle elements, ^^Ne is produced by a- 
captures on to ^''N giving '^^O which subsequently captures 
another a particle to produce ^'^Ne. This takes place in the 
intershell with ^''N being present in the ashes of H-burning. 
Once dredged-up, the ^'^Ne experiences proton-captures in 
the H-burning shell and this afl'ects the abundances of the 
other Ne-Na cycle elements (namely ^°Ne, ^^Ne, ^"^Na as well 
as the two unstable isotopes '^^Na and ^'^Na). The model 
with VW mass loss clearly displays the production of '^^Ne 
from the intershell and the Ne-Na cycle elements from the 
H-burning shell. There is even a trace amount of the radioac- 
tive isotope ^^Na in this model, with a surface abundance of 
4 X 10"^" when the star reaches a temperature of 100,000 K. 

The yields of ^'^Ne and ^^Na have similar uncertainties 
(about 40%). This is because the former comes from the in- 
tershell and all the models dredge-up similar amounts. ^^Na 
is formed from proton captures on to ^^Ne and hence closely 
follows its behaviour. ^^Ne is more sensitive to conditions in 
the burning shell and hence we see a greater uncertainty in 
its yield. In contrast, we see very large uncertainties in the 
abundances of ^^Ne and '^'^Na (by factors of a few), com- 
pared to an uncertainty of 59% in the '^'^Ne abundance. This 
is because ^^Ne is a daughter of the abundant nucleus ^'^Ne 
which does not experience much nucleosynthesis in these 
stars. 

For the heavier elements, the situation is more compli- 
cated. As the star proceeds along the TP- AGB its pulses get 
stronger, with the temperature at the base of the intershell 
convection zone (ICZ) getting higher. This opens up new 
nucleosynthesis pathways. Of particular interest is the reac- 
tion ■^^Ne(a, n)^®Mg, which becomes active at temperatures 
of around 3 x 10* K and provides a source of neutrons. We 
see evidence for this process in the yield and surface abun- 
dance of ^^Mg (and to a lesser extent, ^^Mg). The model 
with Reimers' mass loss has more, hotter pulses and hence 
a greater amount of ^^Mg is produced. The uncertainty in 
the yield due to the choice of mass-loss rate is around 80%. 
This is over twice as great as the uncer tainty associated with 
the ^■^Ne-l-a reaction rates as found bv lKarakas et all (120061 ) 
who studied the production of heavy magnesium isotopes in 
intermediate-mass AGB stars. 



Yields and post- AGB surface abundances for the iron 
group elements are shown in Tables |6] and [T] These elements 
are only affected by neutron captures as their reaction rates 
with charged particles are extremely low at the tempera- 
tures found in AGB stars. The strongest depletion of ^^Fe is 
seen in the VW model as the post- AGB abundances of this 
model are dominated by the last thermal pulse. This is hot 
enough to allow substantial activation of the ^^Ne(a, n)^^Mg 
reaction. The ^''Fe abundance is very well determined be- 
cause this species is very abundant and only a small portion 
of these nuclei suffer neutron captures. The yields are in- 
significant compared to the abundance of ^®Fe because this 
element does not suffer extensive nucleosynthesis on the TP- 
AGB. 

This picture is borne out by the rest of the iron group 
elements. We see evidence for greater processing of these el- 
ements by n-capture in the yield from the Reimers' model. 
This is because this model has the largest number of hot 
(T > 3 X 10* K) pulses. There are large uncertainties in the 
post-AGB abundances of the n-capture progeny of ^^Fe be- 
cause the initial abundances of these elements are low and 
hence small differences in neutron exposure lead to large 
differences in post-AGB abundance. The VW model shows 
the most processing of iron group elements in its post-AGB 
abundances because it is dominated by its last, hot TP. 

The post-AGB surface abundances are more uncertain 
than the yields. This is because the post-AGB surface abun- 
dances depend on what nucleosynthesis is occurring and 
when it takes place relative to the loss of mass. If the enve- 
lope mass is lower when dredge-up occurs, the dredged-up 
material does not get diluted as much and hence the surface 
abundances can change more rapidly. The yields are less sen- 
sitive. The models lose comparable amounts of envelope in 
total so should display similar yields - provided the same 
amount of material is dredged-up to the envelope. 



4 DISCUSSION 

Having discussed the theoretical uncertainties we now move 
on to how the models compare to observations. While we 
could compare our models to observations of AGB star en- 
velopes, we could not be sure how far such a star has gone 
through its TP-AGB evolution. We therefore restrict our 
comparison to planetary nebulae (PNe) and post-AGB ob- 
jects, which we believe have completed their TP-AGB evo- 
lution. 



4.1 Planetary Nebulae 

Some planetary nebulae may be the result of AGB evolu- 
tion, with the material lost during a superwind phase being 
ionised by the hot CO core once it is exposed. Hence PN 
should display the signatures of AGB nucleosynthesis. How- 
ever, we must be careful that the PN we compare come 
from objects of similar mass. For this reason, we do not 
compare these models to nitrogen-rich PN which presum- 
ably come from more massive stars that experience hot 
bottom-burning. iMagrini et al.l (|2004l ) present abundance 
determinations for planetary nebula in the LMC. They find 
e(He) = 10.96, e(N) = 7.46, e(0) = 8.35 and e(S) = 6.81, 



8 R. J. Stancliffe, C. S. Jeffery 



Mass loss 


30Si 


56 pg 


58pe 


eope 


59 Co 


58 Ni 


Reimers 
Blocker 
VW 


1.385(-7) 
1.169(-7) 
1.005(-7) 


-1.604(-7) 
-1.759(-7) 
3.672(-7) 


1.052(-6) 
9.963(-7) 
7.762(-7) 


6.607(-10) 
3.590(-10) 
1.025{-10) 


2.149(-7) 
2.113(-7) 
1.607(-7) 


-1.303(-7) 
-1.221{-7) 
-8.098{-8) 


Difference 
(%) 


3.80(-8) 
32.0 


5.276(-7) 
44.2 


2.76(-7) 
29.3 


5.58(-10) 
149 


5.42(-8) 
27.7 


4.93{-8) 
14.8 



Table 6. Net yield of selected iron group isotoptes in solar masses. For the per cent difference for °°Fe the mean of the absolute values 
has been used as it is more representative. 



Mass loss 


30 Si 


56 Pg 


58pg 


60pg 


59 Co 


58Ni 


Reimers 
Blocker 
VW 


1.020{-5) 
9.785{-6) 
1.363(-5) 


4.563(-4) 
4.625{-4) 
3.882{-4) 


8.058{-6) 
4.883{-6) 
3.919(-5) 


6.439(-9) 
1.929(-9) 
2.647(-8) 


2.688(-6) 
2.047(-6) 
8.813(-6) 


1.858{-5) 
1.916{-5) 
1.211{-5) 


Difference 
(%) 


3.85(-6) 
34.3 


7.43(-5) 
17.1 


3.18(-6) 
197 


2.45(-8) 
211 


6.77(-6) 
150 


7.05(-6) 
42.4 



Table 7. Post-AGB surface abundances by mass fraction of selected iron group isotopes. 



where e(X) — \og{X/ H) + 12 with X and H being abun- 
dance by number. The errors on these measurements are of 
the order of 0.2-0.3 dex. 

We compare the post-AGB abundances of our models to 
these values. The PN that could have formed from our mod- 
els should have come from the final thermal pulses and hence 
have a composition not too dissimilar to the abundances 
of the central object. We therefore have e(He) ~ 11.05, 
e(N) ~ 8.00, £(0) ~ 8.59 and e(S) = 6.88 for the Reimers' 
and Blocker models. Note that the variation between these 
two models is substantially less than the error in the mea- 
surements, being of the order of 0.01 dex. The agreement 
with the observations is reasonable but by no means perfect 
with the nitrogen abundance being somewhat off. This may 
be because we have assume d solar-scaled abu ndances which 
may not apply to the LMC. lRussell fc Dopital (|l992) suggest 
that the nitrogen abundance of the LMC may be a factor 
of 4 lower than the solar-scaled abundance. The VW model, 
which has a significantly lower hydrogen abundance than 
the other two models, gives e(He) = 11.60, e(N) = 8.25, 
e(0) = 9.19 and e(S) = 7.17. Whether these abundances 
would be representative of the PN abundance is question- 
able as there is very little envelope left (after the AFTP) to 
be thrown off from this model. 

Fl uorine has also been recently detected in planetary 
nebula. [Zhang fc Liul l|2005l ') presented measurements of ^''F 
abundances in PN which suggest that in a PN with C/O ~ 2, 
[i9pyi6Qj ^ They con clude that thei r data are consistent 
with the observations of Ijorissen et all l| 19921 ) for giant star 
envelopes. Our models are not consistent these observations 
on two counts. Firstly, our C/O ratios are much higher (by 
about 2.5 in the model with Blocker's mass loss and by about 
4 in the model with Reimers' mass loss). In addition, the 
model [^^F/^^'O] values are too low. The VW model gives 
|i9pyi6Qj _ ]^ g4^ which is in agreement with the measure- 
ments. However, the C/O ratio for this model is an order of 
magnitude too great to match the observations. 



4.2 Post-AGB objects 

We now turn our attention to observations of post-AGB ob- 
jects. Here we must stress that the theoretician's definition 
of a post-AGB object does not match that of an observer. 
To a theoretician, a post-AGB object is an object that was 
recently on the AGB and is now evolving from the AGB tip 
to the white dwarf cooling track. An observer's definition 
of a post-AGB object is a low-mass star that lies above the 
horizontal-branch and to the left of the giant branch. In the 
following, we compare theoretical and observed abundances 
for stars which we believe to be post-AGB stars in the the- 
oretical sense. 

IStasihska et al.l (|2006l ) provide data on C, N, O and S 
abundances for 125 post-AGB objects. We ignore their data 
on R Coronae Borealis (RCrB) stars and extreme helium 
stars as these are unlikely to be post-AGB objects in the 
theoretical sense and hence we do not expect our models to 
fit them. 

Observation of sulphur should be dominated by '^'^S and 
this undergoes little nucleosynthesis on the TP-AGB, being 
subject to limited neutron captures only. Our models suggest 
a value of e(S)= 6.9 — 7.2. At this e(S), the data suggest a 
maximum C/O ratio of about 2.4, whereas our models range 
from 5.1 to 28.0. 

For objects of 0.62 Mq (which is approximately the 
mass of the p o st-AG B objects formed in the models), 
IStasihska et~aLl l|2006t ) find that the total abundance of 
CNO elements over the sulphur abundance (relative to 
solar) should lie between and 0.8. Our models give 
[(C-|-N-|-0)/S]= 0.77 which is consistent with the data. This 
measure should be dominated by the carbon abundance 
in a low-mass object because hot bottom burning, which 
converts carbon to nitrogen, is not active. The models are 
consistent with the oxygen abundance of the observations, 
though there is a large spread observed in the latter. Ob- 
servations suggest £(0) should lie in the range 8.2-9.2. The 
models all fall within this range. The carbon abundances of 
the Reimers' and Blocker models are also consistent with the 



Mass loss and yield uncertainty in AGB stars 9 



data, with e(C) in the range 9.23-9.53. This range is com- 
parable to the data of IStasihska et ahl l|2006l ) who do find 
objects in this range. The VW model is somewhat more C- 
rich. 

We underline that we have only looked at one mass of 
AGB star, whereas the data should contain post-AGB ob- 
jects that have come from AGB stars of a range of masses. 
It is this spread in masses (and other parameters, such as 
metallicity) that should give the observed spread in abun- 
dances. We are also assuming that we would expect to see 
post-AGB objects fro m this mass and meta llicity of AGB 
star in the data set of IStasihska et al] l|2006l ). 



4.3 H-deficient post-AGB objects 

The last thermal pulse of the VW model is qualitatively dif- 
ferent from those of the other two models. Because of its 
low envelope mass, the hydrogen abundance of this model 
falls rapidly during TDUP. While this model does not be- 
come truly H-deficient, it is informative to compare it with 
the H-dc ficio nt post-AGB objects of the PG 1159 class. 
fPreizler fcHebej l|l99^ provide data on the abundance of 
both pulsating and non-pulsating PG 1159 objects. Our 
model is too N-deficient to belong to the pulsating class as 
these objects have N/He ratios of around 0.01. However, it 
has a comparable C/He ratio to the non-pulsators, though it 
is probably too oxygen deficient. The non - pulsa ting PG 1159 
objects in Table 3 of iDreizler fc Hebed (j 19981 ) have C/He 
ratios in the range 0.05-0.3 and O/He ratios in the range 
0.005-0.1. 



4.4 Additional processes 

Whilst the models are consistent with the main features of 
the observations, there are still some discrepancies which 
cannot be explained. The models all produce C/0 ratios 
that are higher than both the observed PN and post-AGB 
object C/0 ratios. There are three possible solutions to this 
problem: 

• TDUP is not as efficient (particularly in the later stages 
of the TP- AGB) as we find, 

• mass loss is stronger in the later stages so that the en- 
velope is removed before the C/O ratio becomes too high, 
or 

• s ome extra mixing process (see e.g. iBoothrovd et al] 
Il995h . which allows material in the envelope to mix to deeper 
in the star, may be active so that partial CN-cycling may 
take place. 

The latter of these options, which would result in a lower- 
ing of the C/0 ratio, is particularly attractive as it would 
help to alleviate the other major problem associated with 
the models, namely the failure to match the observations 
of_[^F/^^0] versus C/0 (see the discussion in lLugaro et al.l 
|2004 for example). 

The failure of the models to match the ^^F observations 
may also be due to the fact that the models do not include a 
^^C pocket. This deprives the models of an important source 
of neutrons and hence may lead to an underproduction of 
fiuorine. 



5 CONCLUSIONS 

We have presented yields and post-AGB surface abundances 
for a 1.5 Mq model of metallicity Z = 0.008, evolved with 
three different mass loss laws. We find the mass-loss law 
has little effect on the mass of the remnant or on the total 
amount of material dredged-up for this mass of star and 
metallicity. We also found efficient third dredge-up even at 
low envelope mass. 

We find a 15% difference in the carbon yields of the 
three models and differences of around 20-80% in the yields 
of many of the other light elements produced in this mass of 
AGB star. Post-AGB surface abundances for these models 
tend to be more uncertain due to dilution effects associated 
with dredging material into envelopes of decreasing mass. 

Comparison of the model abundances to observations 
is mostly favourable, with many of the model abundances 
falling with the observed ranges. However, there are some 
discrepancies. If some additional mixing process were oper- 
ating, partial CN-cycling at the base of the convective enve- 
lope might occur and so could reduce the amount of carbon 
produced by the models. This could go some way to recon- 
ciling the models with the observed C/O ratios, as well as 
observations of ^^F. 

We stress that this work has only examined one mass 
of star and one metallicity and our conclusions only apply 
to these limited circumstances. Future work should expand 
the comparison over a range of masses and metallicities. 
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APPENDIX A: MODEL DETAILS 

Here we present the yields and post-AGB surface abun- 
dances (by mass fraction) of all the isotopes in the nucle- 
osynthesis network. The yields are in Table IXT1 and the post- 
AGB abundances are in Table [K2\ 
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Mass loss 




■*Hc 


12q 


14n 


16 o 


20 Ne 










Reimers 
Blocker 
VW 


-1.380(-2) 

-1.063(-2) 
-2.264(-2) 


2.539(-2) 
2.676(-2) 
1.885(-2) 


4.524(-3) 
4.449(-3) 
3.884(-3) 


3.849(-4) 
3.911(-4) 
3.748(-4) 


1.197(-4) 
1.155(-4) 
1.234(-4) 


1.124(-4) 
1.148(-4) 
1.098(-4) 










Mass loss 




3 He 


■''Li 


''Be 


"B 


13 Q 


14(3 


15n 




18o 


Reimers 
Blocker 
VW 


-1.445(-5) 
-1.475(-5) 
-1.412(-5) 


3.718(-4) 
3.799(-4) 
3.647(-4) 


-2.786(-9) 
-2.844(-9) 
-2.722(-9) 


1.121(-15) 
6.019(-15) 
1.215(-12) 


-1.012(-9) 
-1.033(-9) 
-9.878(-10) 


2.241(-5) 
2.285(-5) 
2.193(-5) 


1.133(-7) 
8.525(-8) 
5.111(-8) 


-6.306(-7) 
-6.423(-7) 
-6.137(-7) 


8.512(-6) 
8.692(-6) 
8.341 (-6) 


-1.441(-6) 
-1.468(-6) 
-1.396(-6) 


Mass loss 


19 p 


21 Ne 


22Ne 




23Na 


24 Mg 


2BMg 


26 Mg 




26A1S 


Reimers 
Blocker 
VW 


3.730(-7) 
3.428(-7) 
2.367(-7) 


3.571(-8) 

3.091{-8) 
1.522(-8) 


2.603(-4) 

2.266(-4) 
1.646(-4) 


0.000 
0.000 
4.283(-13) 


6.301 (-6) 
5.419(-6) 
4.016(-6) 


1.227(-6) 

1.229(-6) 
1.024(-6) 


2.138(-6) 

1.506(-6) 
9.254(-7) 


1.255{-6) 

9.753(-7) 
6.970(-7) 


0.000 
0.000 
0.000 


6.104(-8) 
6.081(-8) 
5.498(-8) 


Mass loss 


^'''Al 


28 g; 


29 gi 


30 Si 




32 g 


33 g 


34 g 


56 Fe 


^'''Fe 


Reimers 
Blocker 
VW 


1.616(-7) 
1.613(-7) 
1.656(-7) 


1.052(-6) 
1.096(-6) 
1.068(-6) 


1.563(-7) 
1.396(-7) 
1.201(-7) 


1.385(-7) 
1.169(-7) 
1.005(-7) 


3.454(-7) 
3.820(-7) 
2.911(-7) 


4.813(-7) 
4.934(-7) 
5.268{-7) 


3.419(-8) 
3.258{-8) 
2.812(-8) 


-2.921(-9) 
-3.249(-9) 
5.055{-9) 


-1.604(-7) 
-1.759(-7) 
3.672(-7) 


5.925(-7) 
5.312(-7) 
4.619(-7) 


Mass loss 


58 Fe 


59 Fe 


60 Fe 


59Co 


68 Ni 


69Ni 


60 Ni 


6lNi 






Reimers 
Blocker 
VW 


1.052(-6) 
9.963(-7) 
7.762(-7) 


0.000 
0.000 
0.000 


6.607(-10) 
3.590(-10) 
1.025(-10) 


2.149(-7) 
2.113(-7) 
1.607(-7) 


-1.303(-7) 
-1.221(-7) 
-8.098(-8) 


2. 745 (-8) 
2.646(-8) 
2.272(-8) 


-2.123(-8) 
-2.224(-8) 
-1.167(-8) 


4.056(-8) 
3.669(-8) 
3.131(-8) 







Table Al. Net yield from the TP- AGB phase in solar masses for all the isotopes in the nucleosynthesis network. All yields are expressed 
as n(m) = n x 10"". 
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^A-Oss loss 




4jjg 


i2r' 


14js[ 


i6n 


20Ne 










Reimers 


6.447(-l) 


3.188(-1) 


2.636(-2) 


9.446(-4) 


4.393(-3) 


7.975(-4) 










Blocker 


6.720(-l) 


3.0.30(-l) 


1.602(-2) 


9.410(-4) 


4.179(-3) 


7.972(-4) 










VW 


3.182(-1) 


5.021(-1) 


1.579(-1) 


8.081(-4) 


7.867(-3) 


7.998(-4) 










Initial abundance 


7.285(-l) 


2.633(-l) 


1.392(-3) 


4.505(-4) 


3.879(-3) 


7.970(-4) 










JVlass loss 


2jj 




7 Li 






i3r< 


14/-' 


15n 




ISO 


Reimers 


0.000 


4.666(-4) 


0.000 


0.000 


4.898(-10) 


4.318(-5) 


2.964(-9) 


8.539(-7) 


1.222(-5) 


6.372(-6) 


Blocker 


0.000 


4.866(-4) 


0.000 


0.000 


5.088(-10) 


4.405(-5) 


1.285(-8) 


8.791{-7) 


1.256{-5) 


6.572(-6) 


VW 


0.000 


2.235(-4) 


6.492(-ll) 


1.193(-9) 


2.387(-10) 


2.717(-5) 


4.613(-6) 


5.339(-7) 


7.701(-6) 


3.827(-6) 


Initial abundance 


1.931(-5) 


1.178(-5) 


3.763(-9) 


0.000 


1.901(-9) 


1.468(-5) 


0.000 


1.755(-6) 


1.563(-6) 


8.718(-6) 


IVlcUSO lOOO 


I9p 


21 Ne 


22 Ne 


22Na 


23 Na 




25 Mp- 


26 Mp- 


26 Al"* 


26A19 


Reimers 


2.859(-6) 


1.947(-6) 


1.884(-3) 


0.000 


5.112(-5) 


2.100{-4) 


4.514{-5) 


4.017{-5) 


0.000 


2.031(-7) 


Blocker 


1.418(-6) 


1.777(-6) 


9.047(-4) 


0.000 


2.955(-5) 


2.082(-4) 


3.384(-5) 


3.471(-5) 


0.000 


1.314(-7) 


VW 


1.436(-5) 


3.128(-6) 


9.985(-3) 


4.215(-10) 


1.772 (-4) 


2.193(-4) 


1.200(-4) 


7.703(-5) 


0.000 


9.671(-7) 


Initial abundance 


1.629(-7) 


1.660(-6) 


5.238(-5) 


0.000 


1.343(-5) 


2.071(-4) 


2.722(-5) 


3.122(-5) 


0.000 


0.000 


JVlass loss 


2^A1 


280: 


29 0; 


300; 


31p 


32q 


33 Q 


34 Q 


56t7„ 




Reimers 


2.348(-5) 


2.619(-4) 


1.448 (-5) 


1.020(-5) 


3.406(-5) 


1.578{-4) 


1.464{-6) 


7.287{-6) 


4.563{-4) 


1.477(-5) 


Blocker 


2.346(-5) 


2.623(-4) 


1.410(-5) 


9.785(-6) 


3.356(-5) 


1.584(-4) 


1.392(-6) 


7.383(-6) 


4. 625 (-4) 


1.326(-5) 


VW 


2.434(-5) 


2.582(-4) 


1.780(-5) 


1.363(-5) 


4.006(-5) 


1.509(-4) 


2.318(-6) 


6.211(-6) 


3.882(-4) 


3.130(-5) 


Initial abundance 


2.332(-5) 


2.627(-4) 


1.378(-5) 


9.463(-6) 


3.281(-5) 


1.592(-4) 


1.296(-6) 


7.507(-6) 


4.703(-4) 


1.148(-5) 


Mass loss 


sspe 


59Fe 


eope 


590o 






eoNi 


61 Ni 






Reimers 


8.058(-6) 


0.000 


6.439(-9) 


2.688{-6) 


1.858(-5) 


1.340{-7) 


7.560{-6) 


5.720{-7) 






Blocker 


4.883(-6) 


0.000 


1.929(-9) 


2.047(-6) 


1.916(-5) 


8.720(-8) 


7.692(-6) 


4.678(-7) 






VW 


3.919(-5) 


0.000 


2.647(-8) 


8.813(-6) 


1.211(-5) 


9.557(-7) 


6.026(-6) 


1.718(-6) 






Initial abundance 


1.487(-6) 


0.000 


0.000 


1.351(-6) 


1.989(-5) 


0.000 


7.877(-6) 


3.457(-7) 







Table A2. Final surface abundances (by mass fraction) of all the isotopes in the nucleosynthesis network. All abundances are expressed 
as n(m) = n x 10"". The initial abundances at the time the star formed are included for comparison. 



